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ABSTRACT 

The broad absorption line (BAL) features seen in a small fraction of quasar optical/UV spectra are 
attributed to bulk outflows away from the quasar core. Observational evidence suggests that dust 
plays a key role in these systems, although whether the inferred dust properties are a signature of 
orientation effects or whether they are indicative of an evolutionary sequence remains an outstanding 
issue. Submillimetre (submm) detections of BAL quasars (BALQSOs), which would clearly help to 
resolve this issue, have so far been sparse. This paper reports on new submm observations of seven 
BALQSOs. The strongest influence on the observed flux is found to be the redshift, with the two 
highest redshift sources appearing intrinsically more submm-luminous than the lower redshift ones. 
Since this trend is also seen in other high redshift AGN, including non-BAL quasars it implies that the 
dust emission properties of these systems are no different from those of the general AGN population, 
which is difficult to reconcile with the evolutionary interpretation of the BAL phenomenon. 
Subject headings: Quasars: Active Galaxies: Ultraluminous Galaxies 



1. INTRODUCTION 

A smaU fraction (~ 0.2 - iBecker et aP 120001: 
iHewett fc Foltzll200,'^ of quasars display broad absorp- 
tion features in their spectra, with widths of several thou- 
sand kms~^, blue- ward of the prominent broad emis- 
sion lines. Seen only in permitted lines, these absorp- 
tion troughs are attributed to bulk outflows from the 
quasar's heart. Typically, the BALs are seen in high 
ionization hues, such as Ly^, Civ A1549 and Ny A1240, 
although more rarely some systems also display low ion- 
ization lines such as Mgn A2800; these latter systems 
comprise ~ 10% of the BALQSO population and are 
referred to as LoBALQSOs, whereas HiBALQSOs dis- 
play solely high ionization absorption features. The BAL 
phenomenon in quasars appears to be part of a much 
broader manifestation of absorption outflows from AGN; 
lower velocity and narrower UV absorption features are 
detected more readily in Seyfert spectra and similar spec- 
tral features from 'associated absorbers' are also detected 
in some (mainly steep-spectrum) rad io-loud quasars (see 
iCrenshaw. Kraemer fc G eorge '20Q2I for a review). In- 
deed, the recently released catalog of BALQSOs from the 
Sloan Digital Sky Survey l)Reichard et aljl2003j) reveals 
numerous quasars possessing weak, difficult to classify 
absorption features, hinting at the possible ubiquitous 
nature of intrinsic absorption outflows. 

The BAL phenomenon, and in particular its relation 
to the overall quasar population, has been the sub- 
ject of debate for a number of years. Two distinct 
competing interpretations for the occurrence of BALs 
amongst quasar spectra are the orientation and evo- 
lution hypot heses. Acc ording to the orientatio n in- 
terpr etation l|Schmidt fc HinesI Il999t IWevmann. et al.l 
ITqqI . all quasars possess BAL outflows with a re- 
stricted covering factor, so that the frequency of de- 
tection simply translates to the rate at which our 
linc-of-sight intercepts the outflow. According to the 
evolution hypothesis ijBriggs. Turnshek. fc Wolfd ITosl 



IVoit. Wevmann. fc Koristalll993(l . the incidence rate of 
the BAL phenomenon is interpreted as the duration of a 
phase of a quasar's natural life cycle. 

Observational evidence supporting the orientation hy- 
pothesis comes largely f rom spectral comparison s of BAL 
and non-BAL quasars IWevmann. et all 119911) and po 
larization studies ( Cohen et al.lll995t iGoodrich fc Millei 



'1995': Hines fc Wills 1995; Hutsemek ers. Lamv. fc Remv 
J,998; Schmidt & Hjncs jL9_99). Evidence in favour of the 
evolution hypot hesis comes largely from IRAS-selected 
BALQSOs (e.g. iBoroson fc Mever^ Il992|) . particularly 
those showing evidence for recent mergers or close inter- 
actions, consistent with proposed evolutionary scenarios 
for quasars (e.g. Sanders ct al. 1988). Further support to 
the evolution hypothesis has been provided indirectly by 
radio observations of BALQSOs, which are inconsis tent 
with simple orientation schemes ijBecker et al.ll2000l) . 

Both hypotheses, however, also rely on observa- 
tions indicating that: 1. BALQSOs ar e substan- 
tially more reddened than no n-BAL quasars ijHall. et alJ 
Il997t iReichard et"an l2003|) : 2. their continua are 
virtually indistinguishable from those of non-BAL 
quasars once dust extinction is taken into account (e.g. 
lYam amoto fc Vanscvicius 1999); a nd 3. LoBALQSOs 
are m uch redder than HiBALQSOs l|Spravberrv fc Foltd 
I1992D . In the orientation scheme, these observa- 
tions support the idea that BAL outflows become in- 
creasingly more dusty as the viewing angle becomes 
more inclined, with LoBALs seen along sight-lines 
grazing the putative nuclear dusty torus which plays 
a key role in orientation-based unification schemes 
for AGN. In evolutionary scenarios, the merging of 
two gas-rich systems fuels both vigorous star for- 
mation and a central AGN, resulting in the strong 
dust emission associated with a n Ultraluminous In- 
frared Galaxy (ULI RG) phase ll^^n^liz^^^tock^gnl 
2f)0Tt ISa,nders et a,].l 17988: Sanders fc Mirabell 1199(1 



Voit. Wevmann. fc Koristall993i) . The overall evolution- 
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SFR 






(mJy) 


(mJy) 


(/i"f?10*MB 


) (/i-^IOI^Lq) (a 


X h7r?Mp,/ 


PSS 1537+1227 


1.20 


-10.1+32.2 


-0.4 +1.0 


<2.1 


<7.2 


<713 


0840+3633 


1.22 


-4.8+9.2 


1.1+1.0 


<2.1 


<7.2 


<713 


1104-0004 


1.35 


7.8+13.2 


0.7+1.8 


<3.9 


<12.8 


<1276 


1556+3517 


1.50 


31.0+28.7 


1.6+1.2 


<1.1 


<3.7 


<371 


1053-0058 


1.55 


6.6+15.3 


0.9+1.2 


<0.7 


<2.5 


<248 


LBQS 0059-2735 


1.59 


19.2+56.3 


10.3+3.3 


7.1+2.4 


23.6+6.9 


2369+764 


Hawaii 167 


2.35 


66.0+20.7 


6.0+1.7 


3.7+1.0 


12.2+3.3 


1219+333 



Table 1 Properties of the sample of BALQSOs, with each column presenting the 
source name, redshift, flux at 450/im and 850/im, the inferred dust mass (md), far 

INFRARED LUMINOSITY (Lfir) AND ST AR FORMATION RATE (SFR), AS DETERMINED VIA THE 

RECIPE of IMcMahon et AL.I (Il999fl . The upper limits represent 3ct values for the 

NON-DETECTIONS . 




h , , , , \ , , , , \ I I I i 

1 1.5 2 2.5 

redshift 



Fig. 1. — The filled circles in this plot represent the bolometric 
luminosities of the BAL systems presented in this study (appropri- 
ate K-corrections have been applied). Error bars represent submm 
measurements only. 



ary path in this scenario is gas-rich mergers — > ULIRG 
LoBALQSO HiBALQSO unobscured quasar. 
Since evolutionary models implicitly predict that 
BALQSOs should exhibit enhanced dust emission with 
respect to non-BAL quasars, one method of distinguish- 
ing between the two competing hypotheses is to measure 
submm luminosities. To date, only a handful of BALQ- 
SOs have been observed at submm wavele ngth s (e.g. 
[Hughes. DunloD, & Rawlings 1997; Lewis e t al.l Il998t 
iPage et al.ll2001l|) . In this paper, we present submm pho- 
tometry observations of seven BALQSOs in order to de- 
termine how their dust emission properties compare to 
that of the general quasar and AGN population. 

2. OBSERVATIONS 

Our target sample was drawn from the literature^ and 
consisted of a range of BAL systems. Given weather con- 
straints, seven of our sample was observed. These were 
PSS 1537+12272, 0840+3633 ( Becker et al. 1997), 1104- 
0004 (Brotherton et al. 199 ^), 1556+3517 JiNaiita et all 
12000^) .' 1053-0058 (Brotherton et all Il998|) and LBQS 

^ This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 

^ Currently unpublished, the position of this quasar was deduced 
from images available in the press as 15 37 42, 12 27 44 (J2000) 



0059-2735 lIHazard et ah l 119871 IWevmann. et a] . 1 11991 . 

whose submm flux are presented in this p aper The de- 
tails o f our ob servations of Hawaii 167 l|Cowie. et alJ 
119941 lEgami e t al. 199^ are presented elsewhere 
(jLewis fc Chapm an 2000). 

This sample was observed in May 2000 with the Sub- 
millimetre Common-User Bolometer Array (SCUBA) on 
the James Clerk Maxwell Telescope on Mauna Kea, 
Hawah. We used the PHOTO METRY three bolome- 
ter chopping mode de scribed in iChanman et all 1)200(1 ) 
and lScott et al.l 1)200(11 to keep the source in a bolome- 
ter throughout the observation. This mode has the ad- 
ditional advantage of allowing a check on the apparent 
detection of a source over three independent bolometers. 
While the 450/im and 850/im arrays are illuminated si- 
multaneously, the bolometer alignment is not perfect, 
and we did not include the 450yLtm offbeams in our final 
flux estimate, except to check that the source had off- 
beam flux consistent with the detection in the primary 
bolometer. 

The observations incorporate chopping (7.8125 Hz) 
and nodding (every 9 seconds), and the final flux den- 
sity in an individual bolometer is a double- difference with 

= 3 beams. The central beam has an efficiency of unity 
and the two off beams have 

. = -0.5exp(-^), (1) 

where d is the angular distance of the off-beam centre 
from the source, and Cb is the Gaussian half-width of the 
beam. For the secondary bolometer the beam efficiency 
is simply 0.5. However, distortion in the field results in 
our chosen third bolometer being slightly offset from the 
source position, resulting in a beam efficiency of 0.44. 
Our 2- and 3-beam measured fluxes agree within < la 
of the primary beam measurement in all non-detection 
cases. For the two sources we claim as detections, H167 
and LBQS0059, the detection signiflcance increases after 
folding in the negative flux density from the two offbeam 
pixels. 

The effective integration time on source varied from 
1200 s to 4800 s for the seven objects in our sample. 

^ The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organ- 
isation for Scientific Research, and the National Research Council 
of Canada. 



Sub-mm emission of BALQSOs 
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Fig. 2. — The filled circles in this plot represent the /gso-optical 
colours for the BAL systems presented in this study. Appropriate 
K-corrections have been applied, so values reflect the rest- frame ra- 
tio, whereas error bars reflect the uncertainty in the submm mea- 
surements only. The open squares in this plot arc a subsample 
of the submm observations of radio galaxies by lArchibald. et ahl 
120011) and are discussed in more detail in the text. 



The secondary was chopped with a throw of 52 arcsec 
to keep the source on bolometer at all times. Pointing 
was checked before and after the observation on blazars 
and a sky-dip was performed to measure the atmospheric 
opacity directly. The rms pointing errors were 1.4 arcsec, 
while the average atmospheric zenith opacities at 450^m 
and 850/zm were 1.7 and 0.21 respectively. The data were 
reduced using the Starlink package SURF (Scuba User 
Reduction Facility, Jcnncss & Liahtfoot 199^ and our 
own reduction routines to implement the three bolometer 
chopping mode. Spikes were first carefully rejected from 
the data, followed by correction for atmospheric opacity 
and sky subtraction using the median of all the array pix- 
els, except for obviously bad pixels and the source pixels. 
The data were then calibrated against standard plane- 
tary and compact Hii region sources, observed during 
the same night. 

3. RESULTS 

Table n summarizes the submm photometry and dust 
properties of our sample of BALQSOs. There is a strong 
redshift dependence upon the submm properties of these 
BAL systems, with only upper limits for all systems 
below z ~ 1.59, despite deep and uniform noise lim- 
its achieved for this sample. The inferred associated 
dust masses, FIR fluxes and star formation rates are 
substantially higher for the two highest- 2; sources com- 
pared to the values calculated from the upper limits 
for the lower-2; sources. Whilst our sample is small, 
these results nevertheless support the notion that the 
high redshift universe is considerably more dusty than 
the local universe. Our results also corroborate other 
separate submm observations of three different BALQ- 
SOs, all with S^mum > 5mJy and a ll lying at z > 1.6 
llHughes. Dunlop. fc Rawlinesl Il997t iLewis et alJ Il998t 
iPage et al.ll27)0ljr 

How do the submm properties of our BALQSO sam- 
ple compare to those of other AGN at high redshift? 
[Priddcv et al. ( 2003) obtained 850/xm fluxes for a sam- 
ple of 57 optically bright quasars in the range 1.5 < z < 
3.0, overlapping with the redshift range of our BALQSO 
sample. With slightly worse noise limits, their study de- 
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tected nine targets at > 3(J significance, a similar detec- 
tion rate to this sample. The BALQSOs presented here 
have 850/im properties very similar to the lPriddev et al~1 
((^03) bright quasar sample, although only a slight trend 
with redshift is found in their sample, and is consistent 
with negative K-correction effects. 

At higher redshifts, however, the dependence of FIR 
luminosity on redshift becomes more evi dent in samples 
of quasars (e.g. Isaak et al. 2002; McMa hon et alJll99flt) 
and e ven more so for radio galaxies (e.g. lArchibald. et alJ 
120011) . In the [Archibald, et al. ( 20Qlil study of 47 radio 
galaxies between 1 < z < 5, 20 of the sources lie in 
the redshift range 1.2 < z < 2.35, overlapping with our 
BALQSO sample, and 4 are detected at 3a, a level simi- 
lar to both our BALQSO sample and the' Priddev et al"! 
(^£)Q^ bright quasar survey. The FIR fluxes of the 
detected radio galaxi es are broadly similar t o those in 
both quasar samples. [Archibald, et al.l l|2001fl . however, 
also found a dramatic increase of FIR flux with redshift 
[^FiR oc (1 -|- z)^], a much stronger evolution than that 
seen in the quasar samples. This trend is interpreted as 
increasingly younger stellar populations associated with 
the radio galaxies at earlier epochs. Fin ally, 850Mm fluxes 
for a sample of X-ray absorbed AGN l|Page et al.ll200Hl 
reveal a very similar trend to that of our BALQSO sam- 
ple, with 4 significant detections above z = 1.7 and only 
upper limits at lower redshifts, with similar detected lev- 
els of 850/j,m fluxes. 

In Figure n bolometric luminosities (total far-infrared 
luminosity plus total UV/optical luminosity) are plotted 
against redshift for our BAL sample. Both the optical 
and submm luminosities have been K-corrected, assum- 
ing fi, oc with a — 0.5 for the optical correction, and 
greybody emission with (3 = 1.5 and T = AOK for the 
submm correction. Infrared luminosity is integrated over 
the greybody from 1100/im to 60/im rest-frame, while 
optical luminosity is integrated over the powerlaw con- 
tinuum from 2OOOA to 5OOOA. The plot also shows just 
the far-infrared luminosities overlaid (smaller symbols), 
showing that the bulk of the bolometric luminosity is 
generated through dust reprocessing, generating ULIRG- 
class luminosities for both submm detections. Note, how- 
ever, that the deep SCUBA limits are insufficient to ex- 
clude the lower redshift BALQSOs from being ULIRG- 
class sources also. The mean submm signal from the 
five non-detections (calculated by weighting each mea- 
surement by its inverse variance) is 0.7±0.5mJy, which 
translates into a far-infrared lumionsity of 8.6x IO^-'^Lq. 
Note also that the negative K-correction in the observed 
850/xm flux results in a quasi-linear flux-luminosity re- 
lation, varying less than 20% over the redshift rang e 
redshifts 1^3 (e.g. Blain. B arnard, fc Chapm an 2003). 
Thus, while the negative K-correction eases the detec- 
tion of distant submm sources, its small variation over 
the redshift range of interest here implies that the de- 
tections/upper limits would have been similar if the red- 
shifts were interchanged. 

Figure El presents a comparison of Ssso/xm-to-R-band 
flux ratios for our BALQSO s ample (filled circl e s) aga inst 
those for a subsample of the [Archibald, et al.l l)2001fl ra- 
dio galaxies for which R-band magnitudes were obtained 
(open squares) . Both the optical and submm fiuxes have 
been K-corrected, and hence represent rest-frame quanti- 
ties. The errors reflect only uncertainties in the submm 
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flux. The strong evolution of submm luminosities ex- 
hibited by these radio galaxies may be a key signature 
of the star formation history of massive elliptical galax- 
ies which are believed to arise from gas-rich mergers. 
Thus, they may offer clues to test the evolution hypothe- 
sis for the BAL phenomenon. The most striking aspect of 
the submm-to-optical flux ratios in Figure [3 is that the 
BALQSOs are substantially less submm luminous than 
the radio galaxies (relative to their R-band flux). This 
is most likely due to increased obscurati on of the radio 
galaxy optical nucleus, since the .Archibald, et all l)200lD 
radio galaxies are steep spectrum sources and thus, ac- 
cording to AGN unification schemes, are probably being 
viewed at large inclinations with respect to their radio 
axes. Note also that that BALQSOs appear to fall into 
a similar distribution as the radio galaxies, with higher 
rcdshift objects possessing more firm submm detections 
than the lower redshift objects, for which only upper lim- 
its could be obtained. Indeed, t his trend prevails in thi s 
entire sample of radio galaxies l|Archibald. et al.ll2001() . 
One final and important point is that there does not 
appear to be a correlation between the submm and op- 
tical fluxes of our BALQSOs (or of the radio galaxies) 
and the high values of Sssofim relative to R-band flux is 
indicative of dust heating by starbursts rather than by 
the BALQSO nucleus. This is also found to be the case 
in several other large quasa r samples (e.g. Ilsaak et all 
120021 iMcMahon et alJll999l and se e also discussion and 
references in lBertoIdifc''Cro3l2flOd) . 

3.1. Selection effects 

It is important to examine whether there are any se- 
lection biases in our BALQSO sample. The BALQ- 
SOs were drawn from quite inhomogeneous surveys; 
e.g. LBQS 005 9-2735 was selected from objective prism 
plates (Hazard et al.lll987tlWevmann. et alill99H) . while 
Hawaii 167 w as found in a spect roscopic followup of K- 
band sources ijCowie. et al.ll993) . Hence, there is no well 
defined selection function. One possible concern is that 
the highest redshift sources are intrinsically the most lu- 
minous, and thus a correlation between total luminosity 
and submm luminosity could be responsible for our bi- 
modal detections between high and low redshifts. How- 
ever, as Fig. 2 demonstrates, there is no evidence for a 
correlation between the submm and optical luminosities 
of our BALQSOs. Although dust obscuration may be 
affecting the submm-to-optical flux ratios, this lack of 
correlation is also readily seen in larg er sam ples of opti- 
cally bright quasars (e.g. Isaak et al] l20fl2l and see also 
references in lBertoldi fc Coxll2002|) . 

4. DISCUSSION & CONCLUSIONS 

This study suggests that the FIR luminosities, inferred 
dust masses and star formation rates of BALQSOs are 
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comparable to those of ULIRG-class sources. Whilst this 
may at first seem consistent with evolutionary models 
for BALQSOs, two key aspects of these submm proper- 
ties, namely the trend of submm flux with redshift and 
the lack of correlation between the submm and optical 
fluxes, are also seen in samples o f non-B AL quasars and 
radio galaxies llArchibald. et al.l f2001: Is aak et al.l 120021: 
IMcMahon et al.lll999|) . If, as is proposed in evolution- 
ary scenarios, BAL outflows are a rapid mass-loss phase 
triggered by a recent gas-rich merger or close interac- 
tion event involving vigourous star formation and asso- 
ciated enhanced dust emission, then it is difficult to un- 
derstand why the submm properties of BALQSOs are 
similar to those of other AGN. While these dust proper- 
ties are consistent with qualitative ideas about the role 
of star formation in the early universe and the evolution 
of the overall AGN population, they do not suggest that 
there is anything remarkable about the dust properties 
of BALQSOs that would be indicative of a connection 
between the BAL phenomenon and the presence of dust. 

If the submm properties of our (albeit small) BALQSO 
sample are verified by larger samples, how do we then 
interpret the BAL phenomenon? It is now clear that 
standard orientation-based models are unable to account 
for all properties of BALQSOs; radio spectral index and 
radio axes measurements indicate that many BALQ- 
SOs simplv cannot be viewed at large inclination angles 
(jBecker et al.ll2000f) . A possible model that can explain 
the radio properties of BALQSOs is one in which the 
BAL material resides in a poorly-coUimated, weak radio 
jet (Kuncic 1999). Such a model can also further explain 
weak, low- velocity absorption outflows in Seyferts, which 
are clearly low-luminosity counterparts to the BAL fea- 
tures in quasars. Again, such outflows are difficult to 
explain with evolution models in which the outflow is 
triggered by a merger; this is simply not the case for 
Seyferts with BAL- like outflows. 

Owing to the small number of BALQSOs in our sam- 
ple, our conclusions are tentative and clearly, larger sam- 
ples are needed in order to verify our results. Ongoing 
intensive surveys, such as the Sloan Digital Sky Survey, 
are proving to be very successful in discovering quasars in 
a relatively unbiased fashion, with the 2dF Quasar Sur- 
vev n ow cataloguing almost 17000 quasars ijBovle et alJ 
|2000|) . providing an ideal basis for such studies. 
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